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An SU{2) KvBLL caloron gas model and confinement 



O 

o 

(N 



cn 
> 

in 
cn 
o 
o 

ix. 



X 



p. Gerhold, E.-M. Ilgenfritz, and M. Muller-Preussker 
Humboldt- Universitat zu Berlin, Institut fiir Physik, 
Newtonstr. 15, D- 12489 Berlin, Germany 
(Dated: July 27, 2006) 

A semi-classical model is developed to describe pure 5(7(2) Yang-Mills gluodynamics at finite 
temperature as a dilute, non-interacting gas of Kraan-van Baal-Lee-Lu calorons including the case 
of non-trivial holonomy. Temperature dependent parameters of the model (asymptotic caloron 
holonomy, caloron density and caloron size distribution) are discussed from the point of view of 
lattice observations and of m-medium modifications of the one-loop caloron amplitude. Space-like 
string tensions running into plateaux at distances R ~ 0.8 — 1.3 fm are obtained and compared to 
lattice results in order to find more precisely the average caloron size. Then, the quark-antiquark 
free energy as predicted by the model is considered. In the confined phase a linear rise with the 
separation can be observed up to _R « 4fm, whereas it runs into plateaux above Tc. Screening effects 
in the adjoint potentials are observed together with an approximate Casimir scaling of the caloron 
contribution to the fundamental and adjoint forces. In Abelian projection, space-like percolation of 
monopoles is found in the confined phase only. Thus, taking the non-trivial holonomy into account, 
confinement properties of pure 517 (2) Yang-Mills gluodynamics can be described by a semi-classical 
approach up to distances one order of magnitude larger than the caloron size. 
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I. INTRODUCTION 

Lattice QCD provides an ab initio description 
of fundamental features of the hadronic world. It 
would be useful, however, to have a continuum 
model of the vacuum structure able to reproduce 
a variety of lattice results and eventually allow- 
ing to compute other quantities not accessible in 
Euclidean lattice QCD. 

What we have in mind here is a reference 
model that could create an ensemble of non- 
Abelian gauge fields, not just various Abelian 
or center projections thereof, analogous to the 
ansatz by Callan, Dashen and Gross U, who 
proposed to describe the QCD vacuum as a su- 
perposition of a dilute set of instantons. This 
approach, widely known as 'dilute gas approx- 
imation' (DGA), was refined by many authors. 
Numerical calculations in such instanton based 
models have been performed by sampling the col- 
lective parameters of their building blocks 
In the case of the 'random instanton liquid 
model' (RILM) ^ the collective parameters were 
sampled directly (independently), whereas im- 
portance sampling was used in the so-called 'in- 
teracting instanton liquid model' (IILM) 0, Q in 
order to account for the residual interactions be- 
tween the pseudoparticles, which are induced by 
temperature and/or the exchange of fermions. 

More recently, superpositions of regular gauge 
instantons, which arc strongly interacting due to 
the slow decay of Af^{x), have been considered for 



which importance sampling is indispensable [g- 
In the present paper we wish to extend the RILM 
in another direction, choosing more general pseu- 
doparticles, calorons with generic holonomy, as 
the basic building blocks. 

Expanding the Euclidean action around such 
a solution and using the method of collective co- 
ordinates 7], it is possible to present the result 
of quantum fluctuations (gauge fixing terms are 
suppressed) as an expression like 



dCJ{C) (det'M(C)) = , (1) 



which is the single-pseudoparticle contribution 
to the partition function. Here det'M(C) de- 
notes the non-zero mode determinant of the 
Hessian M{C) of the Euclidean action S'£;[^], 
parametrized by the collective coordinates (mod- 
uli) C of the classical solution A"^'. Further- 
more, the integrand in Q includes a Jacobian 
J{C), which contains the metric g in the mod- 
uli space in the form y/det g built from the zero 
modes As long as interactions are negligi- 

ble, the integrand in is the probability distri- 
bution used for the direct sampling of collective 
parameters for superpositions of solutions of this 
type. The actual density of such lumps must be 
taken as given by (say, lattice) observations. In- 
deed, a lumpy structure has been observed (and 
its instanton nature taken for granted) in lattice 
studies usin g te chniques like cooling or smooth- 
ing HHIlMllIl or smearing 
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The example of the instanton hquid shows that 
the contributions only of small instantons are 
well under perturbative control by the instanton 
transition amplitude calculated by 't Hooft [T3 |. 
The behavior in the infrared (i. e. the prob- 
ability of large instantons) is described by other 
(mainly classical or Higgs-like) interaction effects 
dealt with in a mean-field fashion or by vari- 
ational techniques _16] . Irrespective of these less 
determined details, the model defines what we 
consider as a model of semi-classical type, in the 
sense that the building blocks of the model are 
classical solutions of the Euclidean equations of 
motion. 

Generically, as the result of mixing of the 
zero modes of individual instantons, the model 
explains the occurrence of a band of near-zero 
modes, i.e. chiral symmetry breaking The 
instanton liquid model is successfully describ- 
ing hadronic correlators and details of hadronic 
structure 0, 0, and is one variant of solv- 
ing the U{1)a problem However, with- 
out adding intricate instanton correlations, the 
instanton liquid model could not describe con- 
finement j22| . This was the motivation to con- 
sider strongly correlated instantons in the regular 
gauge 0. Our paper is aiming to show that the 
lack of confinement can be overcome, for temper- 
atures in the confining phase, by the extension of 
the model from the weakly interacting, singular- 
gauge instanton (caloron) type of classical solu- 
tions to more general solutions with non-trivial 
holonomy [H ll El El . 

This extension of the parameter space of the 
underlying classical solutions can be motivated 
by both: analytical considerations E^, as dis- 
cussed in the following Section, as well as by lat- 
tice observationsjlljlpjpl IM El El El El 
El- In Refs. [20fll3irS'C/(2) lattice gauge 
fields have been analyzed using smearing tech- 
niques and studying their monopole cluster struc- 
ture. Below the critical temperature a part of 
the emerging topological clusters could then be 
characterized to correspond either to non-static 
calorons or to static dyons in the context of 
Kraan-van Baal-Lee-Lu caloron solutions with 
non-trivial holonomy El El El El- The rel- 
ative abundance of topological objects with non- 
trivial holonomy was furthermore shown to be 
temperature dependent and to decrease above 
the critical temperature. In Refs. (3^^] SU{3) 
Monte-Carlo gauge fields have been examined by 
studying the fermionic zero modes. Configura- 
tions with topological charge Q — ±1 probably 
indicating a 3-dyon structure have been iden- 
tified. The 3-dyon structure is a direct fea- 



ture of Kraan-van Baal-Lee-Lu caloron solutions 
in SU{S) with non-trivial holonomy, while it 
can neither be explained as appropriate embed- 
dings of SU{2) instantons into SU{3), nor with 
calorons with trivial holonomy. A full charac- 
terization of SU{3) caloron solutions by cooling 
techniques has been given in Ref. |3| . 

This paper is organized as follows: In Sec- 
tion ^ we will briefly describe the new type of 
classical solutions, calorons of generically non- 
trivial holonomy, to be used in the extension of 
the instanton liquid model. In Section IIIII we 
will describe the problems encountered in the 
construction of superpositions of calorons with 
non-trivial holonomy. In Section IIVI we will 
demonstrate that, under otherwise similar condi- 
tions, non-trivial vs. trivial holonomy determines 
whether the caloron gas confines or not. In Sec- 
tion we will discuss realistic input parameters 
for the model, in particular the analog of the in- 
stanton size distribution, i. e. the distribution 
of the dipole moment of the calorons in terms of 
their constituent monopoles. 

Already here a general remark is in order for 
the orientation of the reader: all this could be re- 
alized in a continuum model. For practical rea- 
sons, however, we will discretize the generated 
gauge field configurations on a suitable grid such 
that the use of lattice techniques becomes pos- 
sible. In Section IVII the simulation results for 
the spatial string tensions and the color averaged 
free energies obtained from our model will be pre- 
sented, before we proceed to Section FVIII where 
we will discuss indirect indicators for confinement 
observed in the monopole structure of the gen- 
erated caloron ensembles (via maximal Abelian 
gauge and Abelian projection). In Section rVIIII 
we shall draw conclusions and give an outlook of 
what should be done next. 



II. CLASSICAL SOLUTIONS OF SU{2) 
YANG-MILLS THEORY 

The instanton, discovered in 1975 is a clas- 
sical solution of the Euclidean Yang-Mills equa- 
tion of motion at zero temperature with localized 
action density, carrying one unit of topological 
charge. For any number of colors, the instanton 
is basically an SU (2) object that is parametrized 
by 4-A^coior = 8 collective coordinates, which are 
its four-dimensional position in space-time, a 
size parameter and three parameters describing a 
global rotation in SU{2) group space. For higher 
A'coior the additional parameters describe the em- 
bedding of SU (2) into SU ( A'^coior) ■ Classical solu- 
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tions with higher topological charge can in prin- 
ciple be constructed by means of the ADHM for- 
malism |37j| . although they are not analytically 
available in general. 

The analog of the instanton in the case of 
non-zero temperature T = 7^ was discovered 
soon after the T — instanton by Harrington 
and Shepard in 1978 "ssj and was accepted as 
an appropriate semi-classical background at fi- 
nite temperature. It will be referred to as "HS 
caloron" throughout the present paper. 

Twenty years later^Kraan and van Baal jlS, 
01 and Lee and Lu [2^ extended the parame- 
ter space of the HS caloron by an additional pa- 
rameter, the asymptotic holonomy Voo which is 
defined as 

7'oo = e'"^^= Jim P(f), (2) 

\x\ — 'OC 

where r denotes the vector of Pauli matrices and 



Fix) = P exp ^ J Ai{x,t)dtj e SU{2) , (3) 

where path ordering is implied by the P sym- 
bol. The generalized solution, which we will call 
KvBLL caloron, is reduced to the HS caloron for 
trivial asymptotic holonomy, taking values in the 
center of the gauge group, Voo G Z{2). In the 
case of general A'coior, this solution can also be 
constructed. It is not a simple embedding of 
an SU{2) solution into the bigger group. The 
present paper, however, is restricted to the case 

of iVcolor = 2. 

Soon after the discovery of the respective solu- 
tions, instanton ensembles as well as HS caloron 
ensembles with trivial holonomy have been used 
to model the QCD vacuum at T=0 or at fi- 
nite temperature, respectively. Concerning non- 
trivial holonomy, the ruling opinion was since the 
beginning of the 1980's that such classical solu- 
tions cannot play a significant role in the QCD 
partition function. This was due to an argu- 
ment given by Gross, Pisarski and Yaffe [s^ that 
gauge fields with non-trivial holonomy are expo- 
nentially suppressed by the perturbative free en- 
ergy j40|. Hence, their potential relevance for 
quark confinement has not been explored so far. 
However, one has to keep in mind that the argu- 
ment does not hold for caloron ensembles with a 
finite density of calorons in space. 

Very recently, Diakonov et al. have calculated 
the holonomy dependence of the free energy of 
a non-interacting KvBLL caloron gas as a con- 
sequence of the caloron quantum weig ht l27^. 



For sufficiently high temperatures, trivial holon- 
omy was shown to minimize the free energy 
and seems therefore dominant in the deconfined 
phase. However, it becomes unstable below a cer- 
tain temperature. Thus, one may conjecture that 
KvBLL calorons with non-trivial holonomy may 
become the relevant degrees of freedom for lower 
temperatures. With this new argument in mind 
and with the KvBLL caloron solution at hand it 
seems a natural step to study KvBLL caloron gas 
models. 

The KvBLL caloron with arbitrary asymp- 
totic holonomy is a self-dual gauge field, i.e. 
F^iv = Ffj^i' with F^^ = ^e^i/a/ji^a/S, solving the Eu- 
clidean equation of motion at finite temperature. 
It was first constructed in the so-called algebraic 
gauge by means of the ADHM formalism 
by Kraan and van Baal and by Lee 

and Lu [2^ as a self-dual gauge field in flat IR'^ 
with the periodicity condition yl°'^(a; -I- f3et) = 
VooA'jj^{x)Vla, where St denotes the unit vector 
along the time direction. The resulting vector 
potential Alf^ defined on can then be trans- 
formed into a self-dual and time-periodic gauge 
field APf^ defined on x K^, by a gauge trans- 
formation ri(x) non-periodic in time 

AP^'^(x,tmod/3) = fl{x)A''Js{x,t)n'<{x) 

- tn{x)d^n^ix) , (4) 

with X ~ (x, t) and 

n{x) = . (5) 

This field is carrying one unit of topo- 
logical charge and the asymptotic holonomy 
Voo- It is the KvBLL caloron in the peri- 
odic gauge with arbitrary asymptotic holonomy 
Voo — exp(27ricJT) S SU{2), and one finds that 
the KvBLL caloron is also described by 8 collec- 
tive coordinates in addition to its holonomy pa- 
rameter uj. However, the physical interpretation 
of these collective coordinates is different from 
that of the instanton parameters, as will be dis- 
cussed in the following. 

One can imagine the KvBLL caloron in the al- 
gebraic gauge as an infinite chain of equally sep- 
arated, identical instantons in flat IR^, aligned 
along the time direction with each subsequent 
instanton rotated by the holonomy Voo relative 
to the preceding one - as long as the size p of 
the instantons is small compared to the inverse 
temperature (3. In that case it is appealing to 
parametrize the caloron in terms of the collec- 
tive coordinates of the instanton. For p <C /3 
the corresponding caloron consists of one lump of 
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FIG. 1: Action density (left) and Polyakov loop distribution (right) of a KvBLL caloron with maximally 
non-trivial holonomy |a3| = 0.25, p = 0.33 fm and 13 = 1 fm, where the action density is given in instanton 
units Si x fm*. Here, pso denotes the 3D-radius of the almost 0(4)-symmetric caloron action lump. 



action, which is approximately 0(4)-rotationally 
symmetric in space-time, and its radius can be 
described by the parameter p. Fig. ^ shows the 
action density and the field of the Polyakov loop 
(-{x), where 
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^color 



Tr P{x) . 



(6) 



This example represents a caloron with max- 
imally non-trivial holonomy and p<ti/3. The 
Polyakov loop goes to zero at spatial infinity ac- 
cording to the maximally non-trivial holonomy. 

In the opposite case, p^P, the caloron is dis- 
sociated into two constituents, as shown in the 
action density plot in Fig. [21 For p — > oo these 
constituents become static in time in an ap- 
propriate gauge, breaking down the prior 0(4) 
symmetry to an 0(3)-rotational symmetry in 3- 
dimensional position space, and the correspond- 

(12) 

ing 3-dimensional radii of the two con- 

stituents converge to constants approximately 
given by p^^^ ~ and Pg'^ ~ 13/uj, respec- 
tively, where w = |w| S [0; 1/2] and uj — 1/2 — uj. 
The constituents can then be identified with BPS 
monopoles [4]| each carrying the fraction 2uj or 
2(D, respectively, of the total action and topolog- 
ical charge. The distance between the monopole 
positions zi, Z2 is given by c?= |zi — Z2I ='Kp'^ / fi. 
Due to the emergence of these two monopoles at 
larger p it seems more natural to use their po- 
sitions zi and Z2 for the parametrization of the 
caloron instead of the instanton parameters in- 
volved in the construction scheme. The corre- 
sponding vector potentials in the algebraic and 



the periodic gauge are then analytically given 
by MM 

+ ^nl.dAncty, (7) 
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+ 4nlud.lncj) + Sf,,i2niUT3, (8) 



with the functions ^{x), ip{x), (^(a;), x(a:), x{x) 
defined according to 

7^2 ^ g2 ^2p4 

'4>{x) = ^ sinh(47rrw) sinh(47rsa;) 

-f cosh{ATTru)) cosh(47rscL') — cos(27rt) 

+ 7rp^s~^ sinh(47rsw) cosh(47rra;) 

-|- 7rp^r~^ sinh(47rr(D) cosh(47rsa') , (9) 

7,2 ^ s2 _ ^2 pi 
ip{x) — sinh(47rra;) sinh(47rsa;) 

-|- cosh(47rr[D) cosh(47rscj) — cos(27rt), (10) 
X(x) = e 



„4^»t^^g-i sinh(47rsa;)e~2^*t 



e*"*"^r-i sinh(47rrw) , (11) 



(12) 
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FIG. 2; Action density (left) and Polyakov loop distribution (right) of a dissociated KvBLL caloron with 
maximally non-trivial holonomy = 0.25, p — 1 fm and (3 = 1 fm, where the action density is given in 
instanton units Si x fm*. For growing p two constituents emerge from the spherical action lump shown in 
Fig. Ifere, pau denotes the radius of the almost 0(3)-symmetric monopole action lumps in 3 dimensions. 



cj^ix) = ^ . (13) 

Here f= x— zi, s — x~ Z2, r — lf], s = \s\, and 
denotes the 't Hooft symbol 14]. For conve- 
nience, the asymptotic holonomy parameter was 
set equal to ujt = lut^, the inverse temperatm'e 
was chosen equal to unity, /3 = 1, and the con- 
stituents were assumed to be separated along the 
z-axis, i.e. z\,Z2 || e-zi with their f-coordinates 
set equal to zero. The general gauge field with 
arbitrary monopole positions, arbitrary tempera- 
ture and holonomy Too — exp(27rja;T) can be ob- 
tained by applying appropriate global transfor- 
mations on the vector potential. These are a 
spatial rotation, a translation in space-time, a 
global gauge rotation O in color space, accom- 
plishing riwTafi^ =a;f , and an adequate rescaling. 
Hence, the KvBLL caloron for a given tempera- 
ture and fixed holonomy can be parametrized by 
the four-dimensional position of its center, the 
size parameter p determining the monopole sep- 
aration, two angles describing their spatial ro- 
tation and one parameter for an residual C/(l) 
gauge rotation around the axis lot in color space. 
This is the parametrization that will be applied 
in the present work. Three further remarks shall 
be given here. 

the case of non-trivial asymptotic holon- 
omy Voo the caloron gauge field does not vanish 
at spatial infinity. In the periodic gauge the A4 
component of the vector potential converges to a 
constant at spatial infinity whereas the remain- 



ing three components go to zero, 

lim AF^^ (x) = 2t:ujt, lim si^) = ■ 

\x\ — *oo — >oc ' ' 

(14) 

The non-vanishing vector potential will become 
an issue when one tries to superpose KvBLL 
calorons. 

(ii) In the case or tD— >0 the asymp- 

totic holonomy becomes trivial, one constituent 
vanishes and the remaining one becomes a HS 
caloron. 

( Hi ) Finally, the anticaloron, which is the anti- 
self-dual analogue of the caloron, carrying nega- 
tive topological charge, can be obtained from the 
caloron gauge field by 

Ar*^(x,i) = -A,{-x,t), i = l,2,3, 
AT'\x,t) = A^{~x,t) . (15) 

III. SUPERPOSITIONS OF CALORONS 

There has been remarkable progress in the con- 
struction of KvBLL calorons with higher topo- 
logical charge [4^ l43l| . Analytical expressions 
for an arbitrary number of calorons with their 
monopoles placed along one axis have been de- 
rived. Analytical parametrizations for the whole 
parameter space of classical solutions with arbi- 
trary topological charge, however, are not avail- 
able in general. Furthermore, configurations of 
mixed self-dual/anti-self-dual character are not 
known as solutions of the equations of motion. 
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just those configurations which are of overwhelm- 
ing importance. Therefore, the model introduced 
in this paper will be based on approximate classi- 
cal gauge field configurations, constructed out of 
single caioron and anticaloron solutions by some 
kind of superposition scheme. In the case of the 
"random instanton liquid" this scheme consisted 
simply of adding N single (anti)instanton gauge 
fields ("sum-ansatz") 

A^(x)=^4)(x), z = l,...,7V, (16) 

i 

where the yljj (x) were chosen in the singular 
gauge, in which the vector potential outside the 
instanton core drops to zero most quickly. Due 
to the non-linearity of in terms of the vec- 
tor potential, this ansatz leads to deviations 
from exact (anti-)self-duality. However, since the 
non-linearity arises from the commutator term 
the superposition approximately de- 
scribes a classical solution as long as at every 
space-time position x the sum of the gauge fields 
is dominated by the contribution A^*'' (x) of one 
single caioron, or if all vector potentials are al- 
most Abelian. 

For instantons, the simple sum-ansatz yields 
good approximations of multi-instanton config- 
urations, unless the separation between the in- 
stantons becomes smaller than the size p of 
their action lumps. The quantitative criterion 
PiPi/ki ~ 2;2|^ < 1/200 for the apphcabihty of 
the sum-ansatz for the superposition of two in- 
stantons with sizes pi, and P2 and positions xi 
and X2 was established in [l5l |. 

When dealing with calorons with non-trivial 
holonomy this criterion no longer holds because 
the KvBLL caioron gauge field does not vanish 
even far away from the monopole locations. At 
first, the potential does not vanish at spa- 
tial infinity in the case of non-trivial holonomy, 
inducing a correlation between calorons even at 
infinite separations. Secondly, there is a string 
("Dirac string") of strong vector potential be- 
tween the monopoles of a caioron with non-trivial 
holonomy. Despite its strength, this string is 
almost free of action, since its vector potential 
is "fine-tuned" . However, there can be strong 
interactions between Dirac strings and further 
calorons. Therefore, more care is required when 
one is going to apply the simple sum-ansatz to 
KvBLL calorons with non-trivial holonomy, in 
order to avoid that the overlapping vector poten- 
tials of two calorons spoil the exact (anti-)self- 
duality of the single solutions. 

The first problem can be overcome quite eas- 
ily by choosing the algebraic gauge for all the 



calorons before adding their potentials. In this 
gauge the potential ^4'^ vanishes at spatial infin- 
ity and the sum-ansatz does no longer yield inter- 
actions between infinitely far separated calorons. 
On the other hand, the vector potentials are then 
no longer periodic. Periodicity has to be restored 
by applying another gauge transformation to the 
sum of the vector potentials. This leads to the 
modified sum-ansatz 

AP^^'^x) = e-2'^^*"^.^AW''^'9(x)-e+2'^**'^^ 

i 

+ 27rtDf • (5^,4. (17) 

However, this superposition scheme is only valid 
under the restriction that only calorons and an- 
ticalorons with identical holonomy be superposed 
because the gauge transformation necessary to 
bring the caioron vector potential from the alge- 
braic to the periodic gauge depends on the holon- 
omy according to l@J. It would not be possi- 
ble to restore the periodicity of a configuration 
of calorons in the algebraic gauge with a global 
gauge transformation, if the configuration were 
consisting of calorons with different holonomies. 

The second problem, involving the Dirac 
strings, will be addressed in an upcoming paper 
and is not discussed here. Although the construc- 
tion, helpful in this case, is of principal inter- 
est, it will turn out in the following Section that 
the calorons considered in this model seem to be 
only weakly dissociated in realistic physical sit- 
uations. They are mostly merged together into 
joined action lumps. This implies that interac- 
tions between Dirac strings and other calorons 
hardly occur. 

As discussed before for the case of instan- 
tons, superposing calorons according to the sum- 
ansatz (|17f) becomes unreasonable if the sepa- 
rations between different calorons or their con- 
stituents, respectively, become too small. In 
such cases an improved superposition scheme is 
needed. A good candidate for such an improved 
scheme will also be discussed in the upcoming pa- 
per. This improved superposition procedure can 
be derived exploiting the ADHM-formalism and 
restores Shuryak's ratio-ansatz |44j when applied 
to instantons. However, for the caioron model, 
introduced in the current work, it will turn out 
that the caioron liquid is sufficiently dilute under 
realistic conditions to justify the sum-ansatz. A 
measure for the quality of a caioron superposition 
is given by the "action surplus" factor 

S^A^ (18) 
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which becomes unity for exactly (anti-)self-dual 
configurations, whereas 7 > 1, if the classical 
equation of motion is piecewise violated. 

So far, the discussion was referring completely 
to the continuum. In order to evaluate a cer- 
tain set of desired observables numerically, it is 
convenient to introduce a grid (called "lattice" 
in the following). In the later calculations, the 
gauge field Af^{x) of the caloron ensemble, which 
is analytically given by the sum-ansatz, will be 
represented by a set of link variables according 
to 

(x+aft \ 
i J A,,{y)dyA . (19) 

The path ordering operator P will approxi- 
mately be accounted for by dividing each path 
of integration into N sub-intervals with their 
lengths a/N<^p being small compared to the 
average caloron size, such that the vector po- 
tential along these sub-intervals can be assumed 
to be constant. In the numerical calculations 
the convergence of the link variables was en- 
sured by adjusting the number of sub-intervals 
dynamically, corresponding to the variability of 
the gauge field, and checked by comparing the 
resulting parallel transporter for different values 
of iV. 



IV. CALORON ENSEMBLES WITH 
TRIVIAL VERSUS NON-TRIVIAL 
HOLONOMY 

Semi-classical studies based on calorons with 
trivial holonomy, i.e. HS calorons, have already 
extensively been undertaken in the past |39l |. 
Apart from Diakonov's new argument 27], open- 
ing the way to consider calorons with non-trivial 
holonomy, i.e. KvBLL calorons, one may won- 
der what the most interesting physical objectives 
are to begin simulations of KvBLL caloron en- 
sembles. A first striking demonstration what the 
consequences of expanding the caloron parameter 
space would be, can be given when we compare 
the heavy quark- antiquark free energy, extracted 
from the Polyakov-loop correlator in a random 
ensemble of HS calorons, with the same quan- 
tity obtained from random ensembles of KvBLL 
calorons with maximally non-trivial holonomy. 

These ensembles have been created on an open 
32'^ X 8 lattice, embedded in a bigger contin- 
uum of the same temporal periodicity that is 
randomly filled with calorons of a certain den- 
sity. The links have been generated as described 



in the previous Section. The inverse tempera- 
ture has been set equal to (3=1 fm. Postpon- 
ing considerations about the realistic choice of 
the model parameters to Section we have here 
assumed a caloron density of n = 1 fm~* and a 
caloron size parameter fixed to p = 0.3S fm. The 
color-averaged heavy quark-antiquark excess free 
energy F°-'"^{R) can be calculated for the gen- 
erated configurations by means of the Polyakov 
loop correlator according to 

F'^-a^R) = -^\n{« £{x) ■ £{y) >v>c) , 

(20) 

where < • >c denotes the average over all config- 
urations and < ■ >v the volume average over all 
pairs of lattice sites x, y with R=\x — y\ fixed. 

Fig.|3jL shows F°''"3 for trivial {u! = 0) and max- 
imally non-trivial holonomy (a; = 0.25), while all 
other parameters are kept fixed. The strik- 
ing observation is that the excess free energy 
of the pair rises monotonously with the dis- 
tance R in the latter case, whereas it runs into 
a plateau, corresponding to the deconfinement 
of quarks, for the HS caloron ensemble in the 
first case. The effect is directly caused by the 
holonomy. This becomes obvious, if one con- 
siders the pair free energy F°'"^ for asymptot- 
ically large distances R. Then the covariance 
covy (f (£"), £{y)) of £(x) and i{y), which is defined 
as <[<e{x)>v - e{x)] ■ [<e{y)>v ~£{y)]>v, goes 
to zero and one obtains 

F-^3(^R)= -1 ln(<covy(£(f),£(y))+ <i>^y>c) 
-i \n{«e>l>c) . (21) 

While < £ >Y= for maximally non-trivial 
holonomy, leading to an infinitely rising quark- 
antiquark potential, the average Polyakov loop 
<£>Y is non-zero for every configuration with 
uj 7^ 0.25 due to the restriction induced by the su- 
perposition scheme, hence forcing the potential 
to run into a plateau. 

The missing .strictly linear rise of the free en- 
ergy on the quark separation even for uj = 0.25, 
which would have been expected at large dis- 
tances R, is an unsatisfactory feature of this spe- 
cial example with fixed caloron size. It can be 
ameliorated by sampling the caloron size param- 
eters p according to a distribution D{p) with a 
finite width. This is demonstrated in Fig. 
where a p-distribution D{p) oc p^/^ exp(— cp^) has 
been used, which will be physically motivated 
later. Here, the constant c is determined by fixing 
the average size parameter p = 0.33 fm. We will 
return to the question of a realistic p-distribution 
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FIG. 3: Heavy quark-antiquark free energy F""®(J?) derived from Polyakov loop correlators measured in 
caloron ensembles (n = 1 fm~*, /J — 1 fm) with trivial (uj — 0) and maximally non-trivial holonomy {lu — 0.25). 
(a) The size parameter is fixed to p = 0.33 fm. (b) The size parameter is sampled according to a distribution 
with finite width. 



in Section IV CI In conclusion of this Section we 
can state that for caloron gases with tempera- 
ture, density and caloron size in the right ball- 
park a reasonable confining potential becomes 
simply a matter of maximally non-trivial holon- 
omy. 



V. THE PARAMETERS OF THE 
CALORON GAS MODEL 

In the following, for a realistic SU (2) KvBLL 
multi-caloron gas model the input parameters 
will be specified to describe pure gluodynamics. 
The model is completely defined by describing 
how the parameters of each single caloron are 
sampled. The 4-dimensional center position, the 
spatial orientation of the two monopoles, a global 
[/(I) rotation around the axis lot in color space 
will be sampled completely randomly, and the 
"size" parameter p will be sampled according to 
a suitable p-distribution D{p,T). The average 
number of calorons placed on the lattice is de- 
termined by the caloron density n(T) and the 
physical four-dimensional volume of the lattice, 
one extension of which is given by the time pe- 
riodicity length (3 = 1/T. Finally, the holonomy 
parameter l3{T), which is assumed to be a func- 
tion of the temperature, determines the type of 
superposed solutions. According to the selected 
superposition scheme introduced in Section Hill it 
has to be the same for each superposed caloron in 
the same configuration. We put uj(T) = uj{T)e3 
without loss of generality. We are aware that 
this setting breaks the global Z{2) invariance. It 
could be easily restored by randomly selecting 



values w(T) or lj{T) = 0.5 - uj{T). 

Expectation values of observables are then ob- 
tained from averaging over the constructed con- 
figurations without additional weighing, since the 
classical action is proportional to the caloron 
density (apart from local violations of the equa- 
tion of motion) which is the same for all con- 
figurations. Quantum fluctuations are approxi- 
mately accounted for in the sampling process by 
an adequate choice of the caloron p-distribution 
Dip,T). 

The physical scale selected by QCD enters the 
calculation through the dimcnsionful parameters 
of the model, which are the caloron density n(T) 
and the p-distribution D{p, T) (both in their de- 
pendence on temperature) . Finally, in the decon- 
fined phase, the temperature dependence of the 
holonomy parameter uj{T) might be important. 
These model parameters should be chosen con- 
sistently with lattice observations. This will be 
described in the following and gives the opportu- 
nity to discuss what is known about these quan- 
tities. We should keep in mind that the confin- 
ing property is independent of the detailed choice 
of parameters, as long as maximally non-trivial 
holonomy is realized in the confined phase. 



A. The holonomy in the confined and 
deconflned phase 

The KvBLL caloron offers the option to set the 
asymptotic holonomy Voo to an arbitrary value 
in SU (2) through its tJ parameter. In a suffi- 
ciently dilute multi-caloron gas the single caloron 
holonomy parameter uj determines the volume- 
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averaged Polyakov loop which is (approximately) 
selfconsistently determining again the holon- 
omy parameter through <£>v = cos{2Tr\uj\). At 
higher densities this simple relation no longer 
holds, and the average Polyakov loop is stronger 
influenced by the internal Polyakov loop profile 
of the solutions. It becomes a function of the 
holonomy parameter uj, the caloron density n, 
and the p-distribution. The procedure in the 
present work is to adopt the holonomy param- 
eter and its dependence on temperature directly 
from lattice results for the renormalized Polyakov 
loop. 

Obviously, the average Polyakov loop on the 
lattice is a bare quantity suffering from ultravio- 
let divergences, since the Polyakov loop is propor- 
tional to the propagator of an infinitely massive 
test quark moving along the time direction. The 
additive mass shift for the test quark resulting 
from the lattice regularization affects the expec- 
tation value of the Polyakov loop by an exponen- 
tial factor 

<\£\>(x exp (-/3to*") = exp (-A^^am*"') (22) 

depending on the length of the loop /? times a 
divergent mass m''™, which is proportional to the 
ultraviolet cutoff 

m''™ oc i , (23) 

where a is the lattice spacing and /3 = A^T-a is 
the inverse temperature. It has been found nu- 
merically in |4It| that the divergent mass term 
^div -g always positive implying that <|^|> de- 
creases with finer lattice discretizations along 
the time direction. For this reason the average 
bare Polyakov loop actually cannot be consid- 
ered as the caloron holonomy parameter. Since 
the Polyakov loop is renormalizable [4^ |43| , it 
is more appropriate to consider a renormalized 
Polyakov loop independent of Nt- 

A possible renormalization procedure for the 
Polyakov loop involving its smtial correlation 
function is presented in Refs. 0,^3' where the 
renormalized Polyakov loop Lj-gj^iT) was defined 
through the free energy of a quark-antiquark pair 
at asymptotic distance. The corresponding re- 
sults are shown in Fig. 0] 

The qualitative statement is that the aver- 
age Polyakov loop is zero in the confined phase, 
refiecting the unbroken global Z{2) symmetry, 
corresponding to a caloron ensemble with max- 
imally non-trivial holonomy, while it gradually 
approaches unity above the critical temperature. 

Finally, Table shows for some selected 
temperatures the renormalized Polyakov loop 
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FIG. 4: The renormalized Polyakov loop in SU{2) 
LGT versus temperature. The Figure is taken from 
Digal et al. ^]. 



T 


Lrcn (T) 


4uj{T) 


<Tc 


0.0 


1.0 


1.10 Tc 


0.58 


0.61 


1.20 Tc 


0.70 


0.51 


1.32 Tc 


0.78 


0.43 


1.54 Tc 


0.85 


0.35 




1.0 


0.00 



TABLE I: Values of the renormalized Polyakov loop 
as obtained by Digal et al. |43 |. The holonomy pa- 
rameter uj{T) for the caloron gas model is fixed (up 
to the Z{2) symmetry between uj and lj) by the dilute 
gas relation Lren(T) = cos{2nu!). 



Lrcn{T) and the corresponding w-parameter in 
the dilute gas approximation. 



B. Choice of the caloron density 

In lattice QCD Monte Carlo studies directly 
examining the topological structure of the glu- 
onic fields one has attempted to get the instan- 
ton (or caloron) density n{T) by counting the 
lumps of topological charge. Since such investi- 
gations in the past required cooling, instantons 
(calorons) with sizes close to the lattice spacing 
have most likely "fallen through the grid" . Fur- 
thermore, the sizes of the observed lumps are lim- 
ited by the lattice size. Therefore, the instanton 
(caloron) densities obtained through this method 
are most likely underestimated. 

An alternative way to determine the caloron 
density uses the topological susceptibility x 
which can be formally defined in Euclidean space- 



10 



time by 



X = J'^^^J d^x<0\f{qix)qm\0> , (24) 



where q{x) is the operator of the topological 
charge density and T the time ordering symbol. 

A first estimate of x in quenched gauge the- 
ory has been obtained from 1 /-/Vroi n r exp ansion. 
The Witten-Veneziano formula [5(1 . l5ll | relates 
the topological susceptibility of pure gluodynam- 
ics to the ry' mass elevation from the masses of 
the pseudoscalar meson octet. This leads to the 
prediction x « (180 MeV)^. 

The connection to the caloron density n is es- 
tablished through the functional integration ap- 
proach, where I24|l can be rewritten as 



d^x q{x) , (25) 



if contact terms are omitted. The total topolog- 
ical charge can be related to the numbers 
of (anti-)calorons in the 4-dimensional volume V 
by Q^N^ — N^, assuming that a collection of 
these objects exhausts the space-time distribu- 
tion of charge. For an uncorrelated caloron gas 
one would expect a Poisson distribution for 
with a variance f7^±=<Af*>, leading to the ob- 
vious result 



n = x 



(26) 



where n = <N'^~\-N^> /V is the caloron- plus- 
anticaloron density. 

From the renormalization properties of the 
Yang-Mills theory, however, it follows that the 
caloron numbers are not Poisson distributed. 
Instead, the dispersion of iV+ + is signif- 
icantly smaller, cr^+_^jY_ = (4/6) • <Af+-|-7V~>, 
where b = llAi"coior/3, resulting from some kind 
of repulsive caloron interaction 0, [E^ . Assum- 
ing that both species, namely calorons and an- 
ticalorons, obey this dispersion relation indepen- 
dently, being still uncorrelated among each other, 
one would get the modified relation (52j | 



n — — y 
4 ^ 



(27) 



between density and susceptibility. This would 
lead to caloron densities approximately twice as 
high for SU{2) than in the naive approach. How- 
ever, despite the existence of caloron interactions, 
we will adopt an uncorrelated sampling of caloron 
positions since the actual correlation between the 



caloron positions remains undetermined. Corre- 
spondingly, we will adopt relation H26|l for our 
uncorrelated caloron gas model to fix the density 
n(T). 

Now the temperature dependence of the topo- 
logical susceptibility shall be discussed. Several 
lattice techniques have been applied to estimate 
X{T). 

The "cooling method" [a 13 evaluates the sus- 
ceptibility H25II by measuring the total charge on 
each independent Monte Carlo configuration by 
means of an adequate lattice operator qL{x), re- 
placing the continuum operator q{x). After some 
cooling (by local minimization of the action) the 
(continuum) topological charge Q of a configura- 
tion can be expressed as Ql— '^j.qL(x) = Z -Q 
with some renormalization constant Z . In the 
confined phase the lattice topological charges Ql 
converge to integer values during the cooling pro- 
cess (Z— s-1), but above the critical temperature 
these plateaux disappear and the method be- 
comes ambiguous. Therefore, this method is not 
adequate to monitor the topological susceptibil- 
ity across the phase transition. 

The "index method" [H HI determines 
the topological charge Q = nL — nR via the num- 
bers of righthanded and lefthanded fermionic 
zero modes. Only Ginsparg- Wilson fermions j56l 
l57j are sufficiently chiral to allow for an unam- 
biguous definition of chiral modes and, conse- 
quently, the topological charge. 

The "field theoretical method" applied in [sl, 
Isgj j evaluates l)24|l by measuring the correlation 
function directly on the lattice, without cooling. 
In this method, the statistical fluctuations could 
be drastically reduced by using improved opera- 
tors for the topological charge density obtained 
from 1 or 2 smearing steps. The lattice suscepti- 
bility XL is connected to the continuum suscep- 
tibility X by 



XL^{Ql/V)=Z^a\ + M 



(28) 



with an additional additive renormalization con- 
stant M arising from the contact terms. 

In Ref. m m the susceptibility x{T) for 
SU{2) has been obtained by the field theoret- 
ical method. Throughout the confined phase, 
at all T<Tc, the susceptibihty x(T) is equal 
to the zero-temperature susceptibility x{T — 0) 
(within error bars). Above the critical tem- 
perature the topological susceptibility drops to 
zero. From the 2-smeared results the value 
x(T = 0) = (198±8MeV)'* is obtained. The 
caloron densities for T>Tc, that will be adopted 
for the later calculations, are based on this mea- 
surement of x{T) and are presented in Table UTI 
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T 


n{T) = x(r) 


<T. 


(198 MeV)* 


1.10 Tc 


(178 MeV)"* 


1.20 Tc 


(174 MeV)"* 


1.32 


(165 MeV)"* 


1.54 Tc 


(157 MeV)"* 


1.79 Tc 


(136 MeV)"* 



TABLE II: Approximate dependence of the caloron 
density n{T) on the temperature T according to 
n{T)=x{T) from I26I I. The topological susceptibili- 
ties x{T) are taken from the SU{2) results obtained 
by AUes et al. 

At this point it should be stressed once more, 
that this estimate of n(T) through x(T) tac- 
itly assumes that there are topological charge- 
carrying objects (calorons) of charge ±1, and no 
correlation exists between the placements of sin- 
gle calorons and anticalorons within the volume 
V . If there were correlations, the simple estimate 
of n{T) would no longer hold. Assuming that 
calorons are exhausting the topological structure 
in the deconfined phase, the decline of the topo- 
logical susceptibility would either imply that the 
caloron density drops to zero or that topologi- 
cally uncharged objects like caloron-anticaloron 
molecules become dominant in the deconfined 
phase, due to an attractive force between oppo- 
sitely charged calorons mediated by the exchange 
of fermions 0, EOf ■ 

Since we are dealing with pure Yang-Mills the- 
ory, no such correlation will be assumed in the 
following, and the decreasing topological suscep- 
tibility x(T') will be plainly interpreted as a de- 
cline of the caloron density n{T). However, re- 
cent lattice observations |33 have given evidence 
that calorons and their monopole constituents 
are not the only type of topological excitations 
in the deconfined phase. Instead, non-self-dual 
magnetic objects seem to contribute to the vac- 
uum structure at high temperatures, above Tc- 
In this respect the model is certainly not com- 
plete to describe the deconfined phase. 



C. Size distributions 

At the classical level the action of every single 
caloron field configuration is So — Sir'^/g'^ inde- 
pendent of its collective coordinates and there- 
fore all calorons would contribute equally to the 
functional integral. However, the classical scale 



invariance is broken at the quantum level by 
the quantum weight as introduced in In 
the caloron model presented here, the quantum 
weight will be approximately accounted for by 
choosing an appropriate size distribution -D(/0, T) 
according to the single caloron quantum weight, 
which is known in 1-loop order. Since the 
problem has already been encountered in semi- 
classical simulations based on instantons, let us 
begin with a short discussion of the instanton size 
distribution. 



1. Instanton size distributions 

The calculation of the 1-loop quantum weight 
in the background of a single instanton in SU{2) 
has been performed by 't Hooft jljj 

Zi"'* = J d^xd^Udpe^^f{p), (29) 

b = yA^color - "^Nf, Co « 0.64191 , 

where the integration J d^x is performed over 
the four-dimensional instanton position in space- 
time and dU denotes the integration over the 
color group SU{2) according to the Haar mea- 
sure. Here, p is the Pauli-Villars mass, i.e. the 
UV cutoff, and g^{fi) denotes the gauge coupling 
given at this cutoff. Since the quantum weight 
is proportional to p''^^, it is divergent in the in- 
frared region for A^coior > 2. Hence, the integra- 
tion over the collective instanton coordinates di- 
verges unless the large p tail is suppressed by 
some other mechanism which would render the 
model finite. 

Several explanations why large instanton 
sizes could be suppressed in the QCD vac- 
uum (confinement effects, higher order effects, 
repulsive instanton interactions) are discussed 
in [gJ. Whereas confinement effects were ex- 
cluded in |63 due to lattice observations, the lat- 
ter two explanations remain viable mechanisms 
for the suppression of large instantons. Assuming 
the running of the coupling constant g becoming 
frozen at sufficiently large p due to higher order 
effects, one obtains a 1/p^ tail [G^I . Such a behav- 
ior was shown to be not in conflict with available 
lattice data. 

A size distribution falling off as 1/p^ would di- 
rectly yield a linear and confining interquark po- 
tential with a string tension proportional to the 



12 



coefRcient in from of the one-over-cube term |63l | , 
but this could not be reconciled with lattice data. 

A third explanation for the large p suppres- 
sion takes the statistical mechanics of an in- 
teracting instanton gas into account. In |15| 
the large distance interactions between instan- 
tons were treated to be of dipole type as pro- 
posed in [J. Furthermore, a hard-core type re- 
pulsion had to be included in order to guarantee 
the consistency of the calculation. This finally 
led to an exponential suppression of large instan- 
ton sizes, motivating the ansatz for the overall 
p-distribution 

D„,st{p)=ap''-''e~'p\ (30) 

The undetermined constants a and c can be fixed 
by the average instanton size and by normalizing 
the p distribution. Exponential suppression fac- 
tors depending on p^ have also been obtained in 
calculations using variational techniques 'l6l and 
in Monte Carlo simulations ,64 | taking only hard 
core type interactions (ssj into account. 

2. Caloron size distributions 

For KvBLL calorons the calculation of the 
quantum weight is more difficult and has been 
performed only recently. Diakonov et al. (2^ 
found an analytical expression for the 1-loop am- 
plitude in the case of an SU{2) caloron. For the 
metric in moduli space see also Ref. |6^. In the 
limit of small p again a factor p^^^ emerges. For 
p3>/3 = 1/T the quantum weight becomes 

= = (31) 

where C« 1.0314. Here, the vectors zi, Z2 
denote the monopole positions, »'i2 = |zi — Z2I 
is the monopole separation and v = At:luT, 
V = AttuT ^ [0; 2ttT] determine the caloron holon- 
omy. P{v) is the 1-loop free energy density for 
a constant gauge field !Q corresponding to u, v. 
As expected, this formula reduces to the quan- 
tum weight of the HS caloron [s^ in the case of 
trivial holonomy. 



For large p, i.e. large ri2, the quantum 
weight is dominated by the exponential expres- 
sion. The polynomial P"{v), appearing in the 
exponent, is of second order and its zeroes are 
v± = 7rT(l ± 1/V3). Therefore, large p param- 
eters are suppressed by the quantum weight as 
long as V <v^ or w > corresponding to nearly 
trivial holonomy, whereas in the opposite case 
the quantum weight diverges exponentially for 
large p. In the latter case a mechanism to cut 
off the p-distribution is needed to guarantee that 
the single caloron remains an important part of 
the vacuum structure. 

For small p, the exponential suppression ansatz 
H30|) reflecting instanton interactions could also 
be applied to calorons. Since this ansatz has been 
deduced for instantons interacting like pointlike 
dipoles, this is a rather crude parametrization 
for the caloron p-distribution because of its more 
complex structure with two emerging monopoles. 
Additionally, the p parameter, which serves as a 
good approximation for the four dimensional ra- 
dius of the instanton action lump, does not de- 
scribe the action lump size in the case of a disso- 
ciated caloron. However, for p<^(3 the caloron 
is not dissociated and the two monopoles are 
merged together in one action lump of nearly 
spherical shape. In this limit p can justifiably be 
interpreted as caloron size. Then it is reasonable 
to take over the distribution (|30|l for sufficiently 
small calorons, provided that the majority of the 
calorons is not dissociated. 

In a lattice study j62| searching for HS calorons 
at finite temperatures, the temperature depen- 
dent average caloron size p(T) has been measured 
by comparing the topological charge density two- 
point correlator with that of a superposition of 
single HS caloron profiles. The conclusion was 
that in the confined phase the average caloron 
size as a function of temperature is approxi- 
mately constant, p(T<Tc) = 0.33 fm, whereas it 
falls off in the deconfined phase. Therefore, the 
same exponential suppression factor exp(— cp^) 
can be applied to calorons at sufficiently small 
temperatures T due to p{T)/f3 ^ 1. 

However, this approximation is only valid up to 
a certain temperature T* < Tc where (3 becomes 
of the same order as p. Above this tempera- 
ture T* the application of the simple instanton- 
radius suppression formula to the (already partly 
dissociated) calorons is no longer justified. In- 
stead, one can consider the statistical mechan- 
ics of an interacting monopole gas in a similar 
manner. Effectively, this can be parametrized 
by the quantity p3D, which is meant to coin- 
cide with p at p-^P (then describing the radius 
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Temperature 


p- distribution 


fixing of coefficients by 


T < Tc 
T > Tc 


D(p,T) = A ■ p"-" ■ eM-cp') 
Dip,T) = A-p>'-'-eM-U^pTf) 


/ D{p,T)dp =1, p fixed 
/ D{p,T)dp =1, p running 



TABLE III: The p-distributions used in the confined and deconfined phase. In the confined phase the coeffi- 
cients A and c are fixed by reproducing a prescribed average caloron size and the normalization of D{p,T). 
In the deconfined phase there is only one coefficient, A, fixed by normalization. 



of the 0(4)-symnietric, non-dissociated caloron), 
and to provide the 3-dimensional radius of the 
caloron constituent monopoles (which have iden- 
tical shape for Ci; = 0.25) at p^p. Connecting 
these two limits continuously, Panip) is a func- 
tion of the caloron parameter p. In terms of this 
quantity the statistical mechanics of an interact- 
ing monopole gas can effectively be described by 
the ansatz 

Di{p, T)dp = DinstiP3D)dp3D (32) 

which should lead to a distribution of the caloron 
"size" parameter p. 

Here, we are using a numerically accessible def- 
inition meeting the above requirements. We de- 
fine P3D as the effective 3-dimensional radius of a 
region where the action density of a caloron (or of 
its monopole constituents, respectively) exceeds 
a certain fraction of the maximal density in the 
caloron (or monopole) center. Its functional de- 
pendence on p is shown in Fig.|31 The precise cut- 
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FIG. 5: Dependence of the effective 3-dimensional 
caloron/monopole radius pso on the caloron param- 
eter p. Here, pso is quantitatively defined in the 
text as the effective, 3-dimensional radius of the re- 
gion where the action density exceeds a certain frac- 
tion of the maximal density in the caloron/monopole 
center. This fraction is chosen such, that pao — p at 

off fraction is chosen such that p^D ~ P for p<^(3. 
Hence, the ansatz (|32|l becomes equivalent to the 
instanton distribution Dinst{p) as long as p<^l3. 
It also remains applicable at high temperatures. 



However, with increasing temperature the extent 
of the monopoles pa^i is kinematically bounded 
by (i/oj. Therefore, this application of l|30|l to 
the interacting monopole gas alone cannot yield 
a cut off for the p-distribution. 

What is needed in order to finally cut-off the 
p-distribution is a suppression of large distances 
d between the constituents. Such an additional 
suppression factor for large "size" parameters 
p 3> /3 can be obtained from the quantum weight 
(|31() . In the 3-dimensional volume V in the 
quantum weight has been interpreted as a specific 
caloron volume (elementary cell per caloron in a 
dense packing) Vcai (p) depending on the size pa- 
rameter p. We introduce a single caloron volume, 
considered in the limiting case of a well dissoci- 
ated caloron and dense packing as the volume of 
a cylinder, circumscribing the caloron 

2 

Vcalip) = Co{io, T)TT0'd, d^^. (33) 

Here (3=1/T sets the scale of the monopole size, 
and d is the monopole separation. Co(w,T') is 
a holonomy dependent, undetermined factor of 
order unity. Substituting the volume V by this 
ansatz for the specific caloron volume Vcai (p) in 
the quantum weight H31I) and neglecting all non- 
exponential terms one arrives at the form 

I?2(p,T)''^°°e-^(^''^)'«(") (34) 

with 

G{uj) = 16Co(w)7r2wV-|-4(w2+w2-4a;w) (35) 

for the large p suppression. For trivial holon- 
omy the polynomial G{uj) becomes unity, and 
the well known temperature suppression for the 
HS calorons 39] is recovered, but for non-trivial 
holonomy the disastrous repulsive (positive) co- 
efficient of d appearing in the exponential is now 
overruled by the specific caloron volume. Choos- 
ing Co(w,T)>0.82 is plausible since Co(w,T) 
was assumed to be of order unity. This is suf- 
ficient to yield an exponentially suppressed p- 
distribution even though in this case no caloron 
interaction (other than the excluded volume in- 
teraction) was taken into account. In principle. 



14 



the coeiBcient Co has to be determined in a self- 
consistent way, such that the sum of all specific 
caloron volumes Vcai equals the total volume V 
up to a factor close to unity. 

The strength of this suppression rises with the 
temperature and depends on the holonomy. Eq. 
(I33|l . approximating the specific caloron volume 
as a cylinder, certainly does not hold for the 
spherical HS calorons well above Tc- There- 
fore, this ansatz is assumed to describe the p- 
distribution at high temperatures but still in the 
confined phase. 

Sewing together both ansatzes corresponding 
to the opposite temperature regimes p<C/? and 
(3 in the simplest form, leads to the ansatz 

D{p,T) = ap'r^'e-'P^o . ^ . ^-iinpTfci^) 

dp 

(36) 

for the /9-distribution in the confined phase. This 
combination of both approaches guarantees that 
the average caloron size is bounded from above 
for all temperatures. To show this, we plot the 
average caloron size p{T <Tc) obtained by this 
ansatz for the plausible scenario u; = 0.25, 
Co = l in Fig. (square symbols). Here c 
was determined such that at zero temperature 
p{T = 0) = 0.33 fm. 

We see, that through this rough ansatz (|36|l 
the average of p at any temperature through- 
out the confined phase, p{T <Tc), is not only 
bounded but also approximately constant as has 
been observed in lattice studies 67]. Here it is 
due to the weakened explicit temperature depen- 
dence for u = 0.25. 

Since the arguments were rather crude and 
the coefficient Co remains undetermined, only 
the ubiquitous form of an exponential suppres- 
sion factor exp(— cp^) will be subscribed for our 
model, while its actual coefficient c should be 
fixed according to lattice observations. In princi- 
ple, these can be obtained by comparing lattice 
measurements of the topological density corre- 
lator with model calculations using the KvBLL 
caloron profile, analogous to the determination of 
p{T <Tc) performed in |67j which were based on 
the HS caloron profile. In the following subsec- 
tion we will find an alternative prescription to fix 
p{T <Tc) based on the space- like string tension 
leading to a self-consistent tuning for our model. 

In the deconfined phase the caloron holon- 
omy ceases to be maximally non-trivial and ap- 
proaches its trivial value as discussed in Sec- 
tion IV Al Therefore, the HS caloron quantum 
weight, which converges to zero for p^oo, ex- 
plicitly determines D{p,T) for T^Tc- Since the 




100 200 300 400 500 



Temperature [MeV] 

FIG. 6: Illustration of the dependence of the average 
caloron size p (solid) and its interval [p — Ap, p — 
Ap] of (standard) deviation (dotted) on temperature 
T. The two different ansatzes for the p-distributions 
for the confined/deconfined phase are continuously 
connected at the critical temperature Tc- For the 
curve indicated by open squares see the text. 

caloron density n{T) is assumed to fall off fast 
beyond the phase transition as argued in Sec- 
tion caloron interactions do not need to be 
considered. Therefore, the coefficient of p^ in the 
exponential is exactly known (set by the temper- 
ature) and the average caloron size arises directly 
from 

Hence, we are left with two different size dis- 
tributions for the confined and the deconfined 
phase, as summarized in Table IIIII For our 
model, we sew together these two ansatzes for 
the p-distribution in a continuous manner, i.e. 

p{Tc)coni. = p{Tc) deconf. 

(37) 

which is shown in Fig. |B1 for the case 
p{T < Tc) = 0.33 fm. We wiU identify the tem- 
perature, at which both ansatzes are continu- 
ously connected, as the critical temperature Tc. 
This establishes a relation between p{T <Tc) 
and Tc that can be used to fix p{T <Tc) in a 
consistent and unambiguous manner within our 
model. 



3. Final determination of p{T <Tc) 

From lattice calculations 69], the zero- 
temperature string tension a{T — 0) in terms of 
the critical temperature is very accurately known 
for SU{2) to be Tc/v/ct(T = 0) = 0.709 (4). Fur- 
thermore, it is well known that the space-like 
string tension crsiT) coincides with the real 
(time-like) string tension at T — 0, but stays 
approximately constant for all temperatures in 
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the confined phase before it even rises quadrati- 
cally with T above the critical temperature [69| . 
Hence, we obtain the relation for the space-like 
string tension 



cTs{T < Tc) 



Tc 
0.71 



(38) 



Together with the relation between Tc and 
p{T<Tc) imposed by H37|l . one obtains a very 
precise determination of p{T <Tc), since the 
space-like string tension <7s{T) is very sensitive 
to the average caloron size parameter p{T). 

Fig. shows the space-like string tension ob- 
tained from our model, as will be described in 
detail in the next Section, for some selected val- 
ues of p{T <Tc) together with those results for 
o's(T'c), that are derived from and JSHJ- 

This approach unambiguously fixes the average 
caloron size to p{T <Tc) — 0.37 fm. It should 
be noticed that this result is not only consistent 
with the lattice observations for the space-like 
string tension ^()9j, but also with the direct lat- 
tice studies for the caloron size l67| . which gave 
p{T<Tc) = 0.33 fm. The smah deviation of ap- 
proximately 10% can easily be explained by the 
fact that the average caloron size in the lattice 
studies was extracted by assuming the underlying 
topological objects to be HS-calorons. It would 
be worthwhile to repeat the calculation of the 
shape of the topological density correlator for un- 
correlated KvBLL calorons and to compare with 
the data. 

As a byproduct of the determination of 
p{T <Tc) we have also fixed the critical tem- 
perature Tc« 178 MeV as well as the zero- 
temperature string tension (t(0)w318 MeV/fm 
within our model. 
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FIG. 7: Space-like string tensions obtained from our 
model for some selected values of p{T<Tc) together 
with the results for as{Tc) derived from 13 711 and 



VI. RESULTS FOR INTERQUARK 
POTENTIALS 

Let us now examine the confinement proper- 
ties of our model at several temperatures. For 
this purpose, caloron ensembles corresponding 
to temperatures T/Tc = 0.8, 0.9, 1.0 for the con- 
fined and r/Tc = 1.10, 1.20, 1.32 for the decon- 
fined phase have been generated and discretized 
on suitable lattices. As in Section IIVI the lat- 
tices do not obey any boundary condition. The 
performed calculations are open volume simula- 
tions, since this is the easiest way to guarantee 
that even strongly dissociated calorons would fit 
on the lattice. The lattice spacing a was chosen 
such, that the ratio p{T)/a stays approximately 
constant for the selected temperatures. The lat- 
tice sizes are listed in Table HVl together with all 
other relevant model parameters, which are cho- 
sen according to Section El Table IIVI also shows 
the measured modulus of the average Polyakov 
loop <K|> together with the input value of the 
average Polyakov loop cos{2ttuj), as well as the 
action surplus factor 7 as defined in lfTH|) . The 
amount of 60 . . . 70% additional action in the con- 
finement phase seems to be acceptable. 

We begin with the discussion of the space-like 
string tension before turning to the more inter- 
esting time-like potentials that will be derived 
from the Polyakov loop correlator. As seen in 
Section IV C 31 the space-like string tension is a 
sensible tool to determine the caloron size pa- 
rameter. 



A. Spatial Wilson loops 

The color singlet ground state energy of a 
heavy quark- antiquark pair can be extracted 
from the Wilson loop 



(39) 



in the limit Rt 00 according to 



E{R) 



lim 



1 



ln<IF(Cfl,fl^) >+Co, 



(40) 

where Cr^r^ denotes a rectangular path with 
spatial and temporal extent R and Rt, respec- 
tively. However, at finite temperature this limit 
can not be realized due to the compactification 
of time. Instead, spafia/ Wilson loops W{Cii^R.^) 
and their associated potential will be considered 
here. It is supposed to coincide with the physical 
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T 


X 


ni [MeV] 




p[fm] 


p/a 


# 


cos{2Truj) 


<\l\> 


7 


0.80 Tc 


32^ X 10 


198 


1.00 


0.37 


2.66 


777 


0.00 


0.13 ±0.01 


1.61 ±0.01 


0.90 Tc 


32^ X 9 


198 


1.00 


0.37 


2.69 


591 


0.00 


0.14 ±0.01 


1.65 ±0.01 


1.00 Tc 


32=* X 8 


198 


1.00 


0.37 


2.67 


526 


0.00 


0.14 ±0.01 


1.69 ±0.01 


1.10 Tc 


32^ X 8 


178 


0.61 


0.33 


2.62 


160 


0.58 


0.43 ±0.01 


1.29 ±0.01 


1.20 Tc 


32^ X 8 


174 


0.51 


0.31 


2.69 


160 


0.70 


0.59 ±0.01 


1.18 ±0.01 


1.32 Tc 


32=^ X 8 


165 


0.43 


0.28 


2.66 


160 


0.78 


0.72 ±0.01 


1.10 ±0.01 



TABLE IV: For the selected temperatures T/Tc = 0.8, 0.9, 1.0, 1.10, 1.20, 1.32 the corresponding model param- 
eters n{T), uj{T), piT), chosen according to Section^ as well as the number of configurations the lattice 
sizes Ns xNt, and the average caloron size in lattice units p/a are listed. Furthermore, the measured Polyakov 
loop <|£|> together with the input value value cos{2nuj) and the action surplus factor 7 as defined in (I18t are 
shown. 



quark antiquark potential at sufficiently low tem- 
peratures due to the effective isotropy of space- 
time at T = 0. For higher temperatures, spatial 
Wilson loops do not provide a physical poten- 
tial. This is obvious, since the spatial string ten- 
sion stays approximately constant in the confined 
phase and even rises beyond the critical temper- 
ature. 

First of all, we check whether the spatial Wil- 
son loops fulfill the area law 



< W{Cr.bJ >cx e 



(41) 



We show in Fig. |H1 the negative logarithm of the 
Wilson loops < W{Cr^r2) > as a function of 
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FIG. 8: Negative logarithm of rectangular Wilson 
loops, —logiW {C R^R^)) J with side lengths R, R2 in 
fundamental (lower curve) and adjoint representa- 
tion versus the area R-R2. The different symbols 
correspond to different side lengths J?2- 

the enclosed area i?-i?2 of the loops. For the 
fundamental representation (lower curve) the de- 
pendence is almost linear as expected, except for 
some loops with strongly deformed geometry, i. e. 
i?2 S> R. This deviation is due to non-area terms 
additionally which have been omitted in the ex- 
ponent of H41|l . For the adjoint representation. 



however, the curve also starts rising linearly at 
small distances but flattens off at larger R. This 
can be seen as a first evidence that screening ef- 
fects for the adjoint charges can be reproduced 
by our model. 

There are various sources for non-area terms 
to appear in the exponent of Eq. (|41|) . In 
general, the effect of various obscuring non-area 
terms can be compensated by extracting the 
string tension from certain products of loops of 
similar shape. A perimeter term, for example, is 
caused by the self-energy contribution, which is 
proportional to the perimeter R+R2, leading to 

< W{Crm2) >« exp(-(Ti?i?2 - a(i? + i?2) - 7)- 
The constant and perimeter terms can be elimi- 
nated by calculating the effective string tension 
from the Creutz ratios iTOll 



a(i?,-R2) - 
1 
2^2 



(42) 



X In 



<VK(Cfi+i.fl,_i)><VF(Cfi_i,fl, + i)> 

<W(Cr.r^)><W(Cr,r^)> 



However, while the extent i?2 of the contour 
Cr^r2 should be sufficiently large to suppress 
contributions of excited states, the expectation 
value of the Wilson loop increasingly suffers from 
statistical fluctuations with growing i?2- Hence, 
a reasonable value for the side length R2 has to 
be chosen to evaluate the effective string tension. 
As a compromise between systematical and sta- 
tistical errors we choose the side length R2 twice 
as long as i?, i.e. i?2 = 2 • i?, which is widely seen 
as a valid, albeit minimal, condition. 

The resulting effective string tensions ob- 
tained from Creutz ratios with i?2 = 2 • i? are 
presented in Fig. for the fundamental rep- 
resentation. For all temperatures shown below 
Tc they run into a plateau of approximately 
296 — 318 MeV/fm within the error bars at mod- 
erate distances of about i?«0.8 fm. For fixed 
p{T) the observed space-like string tension does 
only depend weakly on temperature. An ex- 
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FIG. 9: Effective string tension g(R,Ri) calculated by Creutz ratios (with Ki = 2-R) from spatial Wil- 
son loops versus the distance K at different temperatures r/rc = 0.8, 0.9, 1.0 for the confined and at 
T/Tc^l-lO, 1.20, 1.32 for the deconfined phase in the fundamental (a) and in the adjoint (b) representa- 
tion. 



act temperature independence could easily be 
obtained by a minimal change of the average 
caloron size due to its strong effect on the space- 
like string tension. For the presented tempera- 
tures corresponding to the deconfined phase the 
spatial string tension runs into plateaux at some- 
what smaller distances. The corresponding val- 
ues of the string tension are strongly reduced 
compared to the confined phase. This is due to 
the decreased caloron density, which was chosen 
according to the drop of the topological suscep- 
tibiUty x(2^) the deconfined phase. This is 
a limitation of the current model focussing on 
(anti-)self-dual monopolcs as the exclusive origin 
of the spatial string tension. 

All curves, however, tend to fall off as R in- 
creases. If one takes this noisy observation seri- 
ous it would be interesting to study the contri- 
bution of caloron systems to very large Wilson 
loops analytically. Since the caloron gauge fields 
become Abelian far outside their cores, it is even 
worthwhile and instructive to consider only the 
contributions of the (Abelian) caloron far fields 
of a caloron gas to very large Wilson loops. Such 
analytical calculations are currently under con- 
sideration [Tlj . 

Fig. shows the effective adjoint string ten- 
sions, which are expected to run into plateaux 
before falling off due to adjoint charges being 
screened by gluons. In our model, in the confined 
phase the adjoint string tensions reach plateaux 
of 580 MeV/fm at a distance of i? « 0.8 fm. 
At larger distances the color screening sets in, 
forcing the adjoint force to fall off. We could 
clearly observe the takeover by screening ef- 
fects in the confined phase at distances i? « 1 fm 



when considering "square" Creutz ratios with 
i?2 = R- However, for the extended Creutz ratios 
with i?2 = 2i? the signal becomes very noisy at 
i?> 1 fm, and hence, we can only see slight indi- 
cations for screening effects here. For the decon- 
fined phase the adjoint forces are much smaller 
than below the critical temperature, due to the 
reduced caloron density. 

Finally, the Casimir scaling hypothesis can be 
checked. While the scaling of the forces (ac- 
cording to the quadratic Casimir operators) be- 
tween sources in varies representations of the 
color group is primarily a property of the one- 
giuon exchange, it holds surprisingly well also 
at intermediate distances, thus constraining the 
confining mechanisms. It was an argument 
against the instanton forces that the Casimir 
scaling would be strongly violated (tJ ■ The ratio 
cr''*(i?)/cr-^""''(it!) obtained in our model is plot- 
ted versus the quark separation R in Fig. ^| It 
can be seen that the ratio of adjoint and funda- 
mental effective string tensions begins very close 
to the predicted value at small distances. 

B. Polyakov loops 

With the parameters of the model fixed, we 
are turning now again to the question of confine- 
ment at finite temperature based on the Polyakov 
loop. The color averaged free energy F°-'"^(R) is 
derived from the Polyakov loop correlator by Eq. 
(|2U|) and is presented in Fig. II lb. For the temper- 
atures belonging to the confined phase the short 
distance behavior of the free energy is indepen- 
dent of the temperature. This has also been ob- 
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FIG. 10: Casimir ratio (JAdj/(^Fund of adjoint and 
fundamental forces from space-like Wilson loops 
versus the distance R at the different tempera- 
tures T/Tc = 0.8, 0.9, 1.0 for the confined phase 
and at T/Tc^ 1.10, 1.20, 1.32 for the deconfined 
phase. The theoretical prediction of Casimir scaling 
is marked by the horizontal dashed line. 



served in lattice studies [43 ■ At larger distances 
the curves become linear and dependent on tem- 
perature. What one would expect for the quark- 
antiquark free energy in the confined phase is a 
linear behavior at large distances. 

This is a consequence of a flux tube forming 
between the sources. The flux tube is some- 
times understood as an clastic string. Due to 
higher string excitations of the color flux tube, 
one would expect corrections to this linear behav- 
ior. It remains to be seen whether the caloron gas 
model can generate such a behavior besides the 
linear confinement corresponding to a stiff flux 
tube. 

The next-to-leading order correction 
is a — Ti/f term, the so-called Liischer 
term [t^. Its coefficient is universally given by 
7i = (d — 2)7r/24, where d denotes the space- 
time dimension. This correction is valid for a 
large class of effective bosonic string theories 
at sufficiently large distances. We therefore 
consider the expression 



A{T) 
R 



+ <t{T)- R + C{T) (43) 



with the free parameters A{T), a{T), C{T) as 
the appropriate fitting function for the measured 
free energies in the confined phase at not too 
small distances. Our model is based on long- 
ranged gauge field excitations, and is therefore 
not meant to describe short-ranged physics. As 
a rough estimation we assume the model to be- 
come valid at distances larger than the average 
caloron diameter, i.e. for R>2p{T). We there- 



fore fit the quark-antiquark free energy F"-""^ (R) 
to Eq. (|43|l only for these distances, as shown in 
Fig.lllb. One sees that the measured free energy 
agrees very well with the expected behavior. 

The string tensions obtained from these fits are 
listed in TablelVltogether with lattice results 
for the quark-antiquark free energy obtained for 
SU{2). The comparison with the lattice study 
shows on one hand that our results for the quark- 
antiquark free energy are of the right order of 
magnitude (at least for T/Tc — 0.8, 0.9), but on 
the other hand it also reveals that our model is, 
in its present usage, not capable of reproducing 
the vanishing of the time-like string tension at 
T—^Tc. A possible cure for this limitation of 
the model will be discussed in the conclusions. 

The corresponding coefficients A{T) belonging 
to the 1/r correction term are shown in Table Ivl 
and should be compared with the Liischer co- 
efficient (73j 7L = (d- 2)7r/24 = 0.26. It can 
be seen that the observations are in qualitative 
agreement with the size of the Liischer term. 

While the quark-antiquark free energy rises lin- 
early at large distances for the confined phase, 
the picture changes completely beyond the phase 
transition. Here, the potentials run into plateaux 
providing zero string tension at large distances. 
This effect is due to the holonomy deviating from 
its maximally non-trivial setting. In Section llVI 
we have already seen that the holonomy param- 
eter alone is sufficient to determine whether a 
caloron ensemble provides a confining potential 
or not. The decreasing caloron density certainly 
influences the strength of the obtained potential 
but is not responsible for the plateau building as 
such. For increasing temperatures the holonomy 
parameter becomes closer to trivial and therefore 
the plateau building already sets in at smaller 
distances. 

For the deconfined phase we have again fitted 
the measured potential with H43|l just to check 
the vanishing of the time-like string tension. In- 
deed, almost vanishing string tensions have been 
obtained. Additionally, we have also tried to 
fit the model potentials with some exponentially 
screened expression, since one sees such Yukawa- 
like behavior in Monte-Carlo simulations |4^ . 
but these attempts were unsuccessful. 

Additionally, the adjoint charge- ant icharge 
free energy has been determined and is shown 
in Fig. Illb . Again, the presented potentials co- 
incide at short distances in the confined phase 
and rise approximately linearly up to distances 
i? « 1 fm. At larger distances, the potentials 
run into plateaux due to string breaking as ex- 
pected. In the deconfined phase screening effects 
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FIG. 11: color averaged free energy versus distance R at different temperatures T/Tc = 0.8, 0.9, 1.0 for the 
confined and at T/Tc = 1.10, 1.20, 1.32 for the deconfined phase in the fundamental (a) and in the adjoint (b) 
representation. The fundamental potentials are fitted to 1431 1 for R>2p. This minimal distance is marked. 



T 


a(r)[M^] 


<t(T) 




A{T) 


0.8 Tc 


100.2 ± 2.2 


0.31 


0.44 


0.33 ±0.02 


0.9 Tc 


134.1 ± 1.5 


0.42 


0.29 


0.23 ±0.02 


1.0 Tc 


139.8 ±2.1 


0.44 


0.00 


0.32 ±0.02 



TABLE V: Below Tc'- the string tensions cr{T) and the 1/r fit parameters A(T) as obtained from fitting 
the quark-antiquark free energy to (1431 . cr{T) is presented in physical units and in units of cr(0), where 
(t(0) = 318 Me V/ fm for our model according to Section f V U 31 For comparison, lattice results for a{T)/a{0) 
taken from Ref. |4ai are presented. 
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FIG. 12: Ratio aAdj/apund of adjoint and fun- 
damental forces from Polyakov loop correlators 
versus the distance R at different temperatures 
r/rc = 0.8, 0.9, 1.0 for the confined phase. The the- 
oretical prediction of Casimir scaling is marked by 
the horizontal dashed line. 



are also observed. 

Finally, the Casimir scaling hypothesis shall 
be checked. We shovi^ the Casimir ratio 
cr"*'(i?)/cr (i?) together with its prediction at 
short distances in Fig. E| It can be seen that the 



measured ratio reproduces approximately 85 % of 
the theoretical prediction. Here, the deviation 
from Casimir scaling is much stronger than in the 
previously discussed case of the space-like string 
tension. 



VII. MAGNETIC MONOPOLE 
CURRENTS 

In the dual superconductor picture of confine- 
ment proposed by 't Hooft [7J| and Mandel- 
stam 75], the observation of a linearly rising 
quark-antiquark potential is explained by the for- 
mation of narrow flux tubes of the chromoelectric 
field. Once such field is inserted to the vacuum by 
quark sources it is squeezed under the influence 
of a vacuum "pressure" provided by a Higgs con- 
densate. Furthermore, 't Hooft also conjectured 
that the relevant degrees of freedom responsible 
for the confinement property of this theory are 
actually U{\) degrees of freedom. Therefore, one 
expects the QCD vacuum in the confined phase 
to be populated by random world lines of Abelian 
magnetic monopoles (Dirac monopoles), the ori- 
gin of the Higgs condensate. If there are electric 
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flux tubes, they are enclosed by additional coher- 
ent azimuthal magnetic eddy currents. The nec- 
essary condition for the viability of this descrip- 
tion is the possibility of percolation of Abelian 
magnetic currents. It is therefore of interest 
whether, in the KvBLL caloron gas model, per- 
colation of Abelian magnetic currents can be ob- 
served. 

The Abelian magnetic currents should be ob- 
servable after fixing an Abelian gauge and split- 
ting the gauge field into a (dominating at large 
distances) Abelian component and (important 
only at short distances) non-Abelian compo- 
nents. 

The maximally Abelian gauge (MAG) is the 
most suitable gauge for the purpose of this de- 
composition. In this gauge, even the neglect 
of the non-Abelian degrees of freedom is a rea- 
sonable approximation. The maximally Abelian 
gauge (MAG) is defined as the gauge maximizing 
the gauge functional 

i? = ^rr[r3t/,,^r3!7tj (44) 

by exploiting the non-Abelian gauge freedom. 
Once the extremization is achieved, on average 
all link variables Ux.^i are as close as possible 
to links U'x ^ belonging to the Abelian subgroup 
{exp(i(/)T3) : 0e [0 : 27r]} C SU{2). The projection 

Ux,^l = uq + iuT U'x^^ = exp {i9x,^T3) (45) 

from SU{2) to U{1) is tantamount to sup- 
pressing ui and it2, keeping only the phase 
dx,tj. = arg (uo + ius) as the U{1) gauge field. The 
U{1) gauge freedom of the remaining Abelian 
gauge field Ox,fi, remains unfixed. The Abelian 
monopole currents can be associated with 

the Abelian projected gauge field in the same way 
as in compact QED |7y|. 

Due to magnetic charge conservation 
= 0, the magnetic currents build closed 
loops (clusters) on a finite and periodic lattice. 
In Monte Carlo simulations on sufficiently large 
lattices the network of monopole clusters was 
found to be composed of a single very large 
cluster traversing the whole volume and many 
other clusters of small size. The histogram h{l) 
representing the abundance of connected clusters 
containing monopole world lines of total length 
I (defined by counting the occupied links) was 
found to be split into two very distinctive parts, 
huv{l) for the small clusters and hjii{l) for 
the percolating ones, with a gap separating the 
two Ix^, [ill . The percolating cluster alone was 
shown to reproduce almost the full monopole 



string tension, meaning that the remaining 
small clusters are not relevant for confinement. 
Furthermore, the 3-dimensional (dimensionful) 

density nfj^^ of magnetic currents belonging to 
the percolating cluster has been found to be re- 
lated to the string tension by rif^^ ja^l'^ = 0.65(2) 
for MAG 79]. Although the exact ratio depends 
on the chosen Abelian gauge fixing procedure 
(and is not universal for all possible non-Abelian 
actions) [s^l, the percolation itself of magnetic 
currents, more precisely non- vanishing winding 
in spatial directions, is an indication for quark 
confinement. 

Fig. 1131 shows the average numbers of clusters 
of a certain size per configuration as histograms 
of monopole clusters with respect to their length 
(participating links) in lattice units, h{l), and 
with respect to their maximal 3-dimensional ex- 
tension, iJ(i?(3£')). The histo grams have been 
evaluated for 160 configurations for each of the 
temperatures T/Tc = 0.8, 0.9, 1.0 in the confined 
phase. These measurements have been made on 
periodic lattices, in contrast to the open volume 
simulations described above. Since calorons are 
not periodic regarding the space — like dimen- 
sions, special care has to be taken to minimize 
the resulting action surplus created by bring- 
ing the caloron ensemble onto a periodic lattice. 
We therefore compute for every space-time posi- 
tion X the analytic, single caloron vector poten- 
tial yl^*''(Aa; = x^*^ — x), where x^'^ is the position 
of the i-th caloron, guaranteeing by appropriate 
shift operations that the inevitable discontinuity 
of each caloron's vector potential appears at only 
those space-time positions (on the antipodal side 
of the torus) where the gauge fields are minimal, 
thus minimizing the action surplus. Except for 
this difference, all model parameters correspond- 
ing to these selected temperatures are chosen ac- 
cording to Table Hvl as before. 

Due to the periodicity a subtle definition of the 
cluster extension R^^^^ is needed. At first, the 
smallest cuboid circumscribing the cluster is de- 
termined. If the extent of this cuboid is smaller 
than the lattice size in all directions, then the 
distance between all cluster points is unambigu- 
ously defined. In the opposite case, the cluster 
can be considered as percolating (in one or more 
directions). The observed cluster sizes are much 
larger than the average caloron size p = 0.37 fm. 
A small fraction of the clusters percolate in the 
sense defined above, and some of these clusters 
wrap around the periodic lattice in the space-like 
directions, carrying a non-zero winding number. 
As discussed above this is a signal for confine- 
ment. 



21 




0.01 H 



0.001 



10 20 30 40 50 60 70 80 90 100 
Cluster length in lattice units 




0.01 H 



0.001 



10 21) 10 .50 60 TO 8(1 90 100 

Cluster length in lattice units 




0.01 



0.001 



0.01 



0.001 



0.01 



0.001 



0.01 



0.001 



10 20 30 40 50 60 70 80 90 100 
Cluster length in lattice units 




2 3 4 5 6 
3D-extent of cluster [fni] 



FIG. 13: Histograms showing the average number of clusters per configuration of a certain cluster size, in 
(a) characterized by the number of occupied links in the cluster, in (b) characterized by the 3D extent of the 
cluster, for the temperatures T/Tc^ 0.8, 0.9, 1.0 in the confined phase. Link numbers larger than 100 are set 
equal to 100. 



In contrast to that, the pictures looks com- 
pletely different for temperatures above Tc- 
Fig. El shows the histograms of monopole clus- 
ters with respect to their length in lattice units, 
and with respect to their maximal 3- 
dimensional extension, H{R^^^^). Again, 160 
configurations have been evaluated for each of the 
temperatures T/Tc = 1.10, 1.20, 1.32 in the de- 
confined phase. One sees that both distributions 
are strongly contracted towards smaller cluster 
sizes. No percolation or wrapping around the 
lattice is encountered anymore for temperatures 
above the critical temperature. This observation 
is in agreement with the breakdown of (time- like) 



confinement that we had found from the study of 
the Polyakov correlators. 

However, it should be noted that the actual 
sizes / of the percolating clusters are one or two 
orders of magnitude smaller than those seen in 
the Monte Carlo ensembles |83|. This is only 
natural since the presently evaluated configu- 
rations are constructed by superposing single 
calorons. Each non-dissociated caloron gives rise 
to a closed cluster traversing through both of its 
constituents such that this cluster entirely lies 
within the caloron core. For dissociated calorons 
two disconnected monopole clusters (world lines) 
emerge wrapping around the periodic lattice in 
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FIG. 14: Histograms showing the average number of clusters per configuration of a certain cluster size, in 
(a) characterized by the number of occupied links in the cluster, in (b) characterized by the 3D extent of the 
cluster, for the temperatures T/Tc = 1-10, 1.20, 1.32 in the deconfined phase. 



time direction. Since the caloron gauge field is 
non-Abelian only in its core, and the vector po- 
tential away from the core can be approximately 
identified by the Abelian projection, the result- 
ing multi-caloron gauge field is Abelian almost 
everywhere. Since the calorons are sitting rel- 
atively close, the monopole world lines reorga- 
nize themselves into a percolating network with- 
out essentially increasing the density of monopole 
currents. Hence, the observed monopole cluster 
sizes in the confinement phase are barely exceed- 
ing the minimal cluster sizes necessary to facili- 
tate the required monopole percolation. This can 
be seen, if one compares the average monopole 
cluster length in lattice units, as measured in 



the simulation, with the average cluster length, 
that would be obtained, if each single caloron's 
monopole cluster was not effected by the presence 
of the other calorons. Table IvTI shows simulation 
results together with the estimates for the av- 
erage cluster length, as obtained from the given 
p-distribution with p — 0.37 fm. It can be seen 
that the total length of the monopole clusters 
induced by the single calorons is increased only 
by 25 — 35% in order to unite closely situated 
monopole clusters to form a connected monopole 
cluster. 
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T 


Measured 


Estimated 


0.8 Tc 


14.7 


11.5 


0.9 Tc 


14.1 


10.9 


1.0 Tc 


12.4 


9.2 



TABLE VI: Average cluster lengths in lattice units as 
measured in the simulations. The result is compared 
to an estimate given by the caloron p distribution, 
assuming each single caloron monopole cluster to be 
unaltered by the presence of other calorons. 



VIII. SUMMARY AND CONCLUSIONS 

We have introduced a semi-classical model 
meant to describe the Euclidean field histories 
of (quenched) QCD as a dilute gas of KvBLL 
calorons with adjustable holonomy. We have 
shown, that this has the striking impact that 
the confinement/deconfinement phase transition 
can in principle be modelled by varying noth- 
ing but the holonomy parameter, starting from 
its maximally non-trivial value ulT <Tc) = 0-25 
for the confined phase, towards its trivial setting 
u){T^Tc)^0 or 0.5, respectively, as the temper- 
ature increases, such that very high above the 
critical temperature the well-known HS-caloron 
is finally recovered. 

In order to work out the model more quantita- 
tively and in detail, the parameters of this model 
have been chosen in agreement with lattice obser- 
vations. Firstly, the holonomy parameter u{T) 
was fixed in accordance with lattice studies of the 
renormalized Polyakov loop in SU{2) by exploit- 
ing the low-density relation <|£|>=cos(27r|a;|). 
The reliability of the current determination of 
uj{T) could be improved, if it would be deter- 
mined in a self-consistent manner. However, as 
long as the holonomy is maximally non-trivial 
in the confined phase and becomes trivial only 
at high temperatures above Tc, the exact func- 
tional form of w(T) has only little infiuence on 
the observed "asymptotic" string tensions. 

Second, the caloron density n{T) has been 
determined by lattice results for the topolog- 
ical susceptibility x(^) by assuming calorons 
and anticalorons to be uncorrelated. Since we 
do not take any caloron/ anticaloron correlation 
into account for the sampling of their positions, 
this is the adequate assumption with respect to 
our model, although calorons and anticalorons 
are actually known to be correlated. This has 
lead us to adopt a constant caloron density 
n(T) = l fm~^ as usual in the confined phase 
and a decreasing density oc x(T) in the decon- 
fined phase. This is not compelling since the 



decreasing topological susceptibility could also 

be explained by topologically uncharged objects 
becoming dominant, such as caloron-anticaloron 
molecules, although there is no independent in- 
dication for that to happen in quenched QCD. 

Also in another respect the current model is in- 
complete in the deconfined phase because it does 
not consider non-perturbative configurations of 
other origin. This incompleteness becomes ob- 
vious by seeing the space-like string tension de- 
creasing above the critical temperature due to 
the reduced caloron density. This is unphysical 
since the space-like string tension is known to 
rise at higher temperatures, completely indepen- 
dent from the topological susceptibility which be- 
comes tiny. Apart from the loophole forcing us to 
accept a decreasing density of calorons, it could 
well be that monopoles being neither self-dual 
nor anti-self-dual build the spatial string tension 
throughout the deconfined phase. 

The breakdown of the physically confining po- 
tentials above the critical temperature in our 
model is driven by the holonomy alone. The de- 
crease of the caloron density would not be re- 
quired. 

Third, the /9-distribution was determined by 
reinterpreting the single caloron quantum weight 
in the dense-packing limit. Due to its diver- 
gence at p^oo for non-trivial holonomy, ad- 
ditional arguments rendering the p-distribution 
finite had to be invoked. Since its analytical 
form could only be qualitatively motivated, a 
more sound investigation of the p-distribution 
in a dense medium of calorons would be worth- 
while. However, having fixed the analytical form 
of the size distribution, the average caloron size 
in the confined phase has been unambiguously 
determined by comparing the space-like string 
tensions with those obtained in lattice studies. 
Hence, the usual criticism of instanton/caloron 
based models, namely that almost every desired 
result could be obtained from this kind of mod- 
els by suitably setting the average caloron size, 
is actually not applicable here. 

The color-averaged quark-antiquark free en- 
ergy was then extracted from the Polyakov loop 
correlation in multi-caloron configurations con- 
structed according to a simple superposition 
scheme. Since our model is based on long-ranged 
gauge field excitations of average diameter 2p(T), 
the model is not supposed to yield reasonable re- 
sults at distances smaller than R<2p{T). For 
R>2p{T) the obtained free energies show the 
expected linear plus Liischer term behavior up 
to distances of i?w4 fm in the confined phase. 
Above the critical temperature the obtained po- 
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tentials run into plateaux corresponding to the 
deconfinement of quarks. The obtained string 
tensions are of the correct order of magnitude as 
can be seen by a comparison with lattice stud- 
ies. However, our model is, so far, not capable 
to reproduce the correct temperature dependence 
of the time-like (color-averaged) string tension, 
which is known to gradually decrease with ris- 
ing temperature before it is going to vanish at 
T = Tc- This limitation is caused by the fact that 
our model does not allow for mixed holonomies 
within a single gauge field configuration. A mini- 
mal improvement of the model would be to sam- 
ple the uniform caloron holonomies in the con- 
fined phase according to a distribution with non- 
vanishing, temperature dependent width (in ac- 
cordance to what is seen in lattice Monte Carlo 
simulations) rather than being fixed exactly to 
zero [uo — 0.25). 

By studying the potentials in the adjoint rep- 
resentation screening effects due to string break- 
ing could be observed. For small distances a ra- 
tio between the fundamental and adjoint effective 
string tensions of cr^'^-'/cr-^""''' ^ 2.2 was obtained, 
which is close to but smaller than the ratio 8/3 
predicted by the Casimir scaling hypothesis. For 
space-like string tensions screening effects could 
also be seen and the Casimir hypothesis was even 
better fulfilled. Here, a ratio cr'^^/cr-'^""'' « 2.6 
could be observed for all selected confining tem- 
peratures. 

Finally, the formation of large magnetic 
monopole clusters was studied in the caloron gas 
model. Percolating clusters traversing the whole 
lattice volume were observed for the confining 
phase, which is a necessary condition for con- 
finement. For the deconfined phase the distribu- 
tion of monopole cluster sizes is contracted to- 
wards smaller cluster sizes and no spatially per- 



colating clusters are found anymore, which is in 
agreement with Monte Carlo lattice observations. 
Not unexpectedly, the observed monopole clus- 
ters are one or two orders of magnitude smaller 
(in path length) than those obtained in Monte 
Carlo studies. The average cluster lengths in the 
caloron ensembles are only slightly increased by 
rearrangements induced by the caloron superpo- 
sition compared to the sum of the cluster sizes of 
the separate calorons. 

In conclusion we would like to say, despite 
the limitations that our model still suffers from, 
the KvBLL caloron gas model offers an interest- 
ing scenario for the confinement/deconfinement 
phase transition. With the proposed improve- 
ments, the KvBLL caloron gas model most prob- 
ably will give a consistent picture of the phase 
transition. Therefore it is a worthwhile topic to 
further investigation. 
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